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UNSWEPT  WINGS  AT  ZERO  LIFT 
By  Jack  N.  Nielsen 


SUMMARY 

The  linear  theory  for  determining  the  pressure  distribution  at 
supersonic  speeds  on  wings  of  symmetrical  section  at  zero  lift  has 
been  used  to  calculate  the  pressure  drag  coefficients  at  zero  lift 
of  the  family  of  unswept,  untwisted  wings  with  the  diamond  profile. 

On  the  basis  of  the  method  of  R.  T.  Jones,  which  was  presented  in 
NACA  TN  No.  1107,  a  general  expression  has  been  found  for  the  drag 
coefficient  at  any  supersonic  Mach  number,  aspect  ratio,  and  taper. 
The  general  expression,  which  is  too  unwieldy  for  use  in  analysis, 
has  been  used  in  the  construction  of  nondimens ional  charts  permitting 
a  rapid  estimation  of  the  drag  coefficient  for  a  given  wing  at  any 
Mach  number.  For  wings  of  diamond  or  rectangular  plan  form  the 
general  expression  reduces  to  simple  formulas  for  the  drag  coeffi¬ 
cients.  The  nondimensional  charts  indicate  that,  at  very  low 
aspect  ratios,  wings  of  diamond  plan  form  have  the  least  drag  and 
that  increasing  the  ratio  of  tip  chord  to  root  chord  increases 
the  drag  markedly.  However,  at  large  aspect  ratios  rectangular 
wings  have  the  least  drag  and  decreasing  the  ratio  of  tip  chord 
to  root  chord  increases  the  drag  slightly. 


INTRODUCTION 

Methods  for  determining  the  aerodynamic  characteristics  of 
various  aircraft  components  at  supersonic  speeds  have  recently 
become  the  subject  of  widespread  attention.  Of  particular  interest 
are  methods  which  permit  determination  of  wing  characteristics  in 
three-dimensional  flow.  Substantial  progress  in  the  development 
of  a  theory  for  the  three-dimensional  supersonic  wing  has  been 
made  by  linearizing  the  problem  on  the  basis  of  small  perturbation 
theory.  This  makes  possible  the  use  of  the  superposition  principle 
and  enables  separation  of  the  pressures  acting  on  an  uncambered 
wing  into  (l)  a  pressure  distribution  due  to  thickness  which 
occurs  at  zero  angle  of  attack  plus  (2)  a  pressure  distribtuion  due 
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to  angle  of  attack  which  occurs  on  a  flat  plat©  of  the  wing  plan 
form  at  the  wing  angle  of  attack.  The  pressure  distribution  due 
to  thickness  for  thin  symmetrical  wings  the  surfaces  of  which  are 
formed  by  planes  can  be  determined  by  the  method  of  R.  T.  Jones 
(reference  l) .  Essentially  the  method  consists  of  orienting 
oblique  line  pressure  sources  and  oblique  line  pressure  sinks  of 
different  strengths  in  a  manner  to  satisfy  the  boundary  conditions 
and  then  adding  their  pressure  fields.  An  equivalent  method  has 
been  published  by  Puckett  (reference  2) . 

The  foregoing  methods  have  been  rather  extensively  applied  to 
wings  of  triangular  plan  form,  but  as  yet  only  limited  study  of 
the  effects  of  changing  plan  form  by  varying  aspect  ratio,  taper, 
and  sweep  has  been  carried  out.  (Reference  3  is  suoh  a  study.) 
Accordingly  in  this  report  the  method  of  Jones  has  been  applied  to 
determining  the  pressure  drag  coefficient  at  zero  lift  of  a  family 
of  untwisted  wings  of  diamond  profile  differing  in  aspect  ratio  and 
taper.  For  all  wings  the  line  of  maximum  thickness,  that  is,  the 
mid-chord  line,  is  unswept.  Tapers  from  zero  to  unity  have  been 
considered  for  all  aspect  ratios. 


SYMBOLS 


A 

8 

C 

CD 

M 

P 

Po 

P 

*o 

S 

t 

t 

c 


aspect  ratio 
semispan,  feet 


root  chord,  feet 

wing  pressure-drag  coefficient 

free— stream  Mach  number 

local  static  pressure,  pounds  per  square  foot 


free _ stream  static  pressure,  pounds  per  square  foot 


pressure  coefficient 


f re ©-stream  dynamic  pressure. 


pounds  per  square  foot 


wing  area,  square  feet 

wc-Hrmun  thickness  of  wing  section  at  root  chord,  feet 

wing  thickness  ratio,  fraction  of  chord,  approximately  equal 
to  tangent  of  semivertex  angle  of  wing  section  measured  in 
streams! se  direction 
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A  sweep  angle  of  leading  edge,  radians 

6  polar  angle  measured  from  downstream  direction,  radians 

a  tan  0 

t  tan  A 

E.P.  real  part  of  a  complex  quantity 
X  ving  taper,  ratio  of  tip  chord  to  root  chord 

Subscripts 

a  leading-edge  sources 

b  midchord— line  sinks 

c  trai ling-edge  sources 

d  leading-edge  image  sinks 

e  midchord— line  image  sources 


ANALYSIS 
Types  of  Wings 

For  the  purpose  of  the  analysis  it  has  been  convenient  to 
subdivide  the  wings  under  consideration  into  eight  types  as  shown 
in  figure  1.  From  this  figure  it  can  be  seen  that  certain  relation¬ 
ships  between  aspect  ratio  and  taper  define  boundaries  which  determine 
the  wing  type.  Along  boundary  III  the  leading  and  trailing  edges  are 
coincident  with  their  respective  Mach  lines.  Above  this  boundary  the 
leading  and  trailing  edges  are  swept  in  front  of  the  Mach  lines 
emanating  from  their  foremost  points  and  are  termed  supersonic  leading 
and  trailing  edges.  Below  this  boundary  the  leading  and  trailing 
edges  are  swept  behind  their  respective  Mach  lines  and  are  termed 
subsonic  leading  and  trailing  edges.  This  figure  may  be  used  for  any 

Mach  number  by  using  kjw— 1,  the  effective  aspect  ratio,  in  place 
of  the  aspect  ratio. 

In  wings  of  type  1,  the  Mach  lines  all  intersect  the  trailing 
edges  of  the  adjacent  half  wing.  In  wings  of  type  2,  the  Mach  lines 
from  the  extremities  of  the  leading  edges  intersect  the  trailing 
edges  of  the  opposite  half  wing,  and  the  Mach  lines  from  the  leading 
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vertex  intersect  the  wing  tips  at  a  point  between  the  midchord  line 
and  the  trailing  edge.  The  differences  between  wings  of  other  types 
can  easily  be  seen  in  figure  1. 


Determination  of  Wing  Pressure  Field 

The  wing  pressure  field  can  be  built  up  by  superimposing  a  number 
of  simple  pressure  fields  which  are  basically  of  two  types.  These 
types  of  pressure  fields  are  those  given  by  Jones  in  r®*®*®11®®  }  f?F  . 
a  semi-infinite  symmetrical  wedge  swept  in  front  of  or  behind  the  Mach 

line  from  its  apex. 

For  a  semi-infinite  symmetrical  wedge  with  a  supersonic  leading 
edge  (fig.  2(a))  the  pressure  field  is  conical  and  for  Vr  =  2  is 
given  by  the  following  equation: 


E.P. 


2  t/c 


cos 


— l 


tan  A  -  tan  9 


'  —  —  .  —  WW 

„t/l  _  tacrA  |l  -  tan  A  tan  e 


(1) 


Between  the  leading  edge  where  9  +  A  * 
tan  0=1  this  equation  reduces  to 

2  t/c 


2L  and  the  Mach  line  where 
2 


4  1  —  tan2  A 

which  is  the  Acker et  equation,  as  given  by  reference  4,  for  the 
pressure  coefficient  based  on  the  Mach  number  component  perpendicular 
to  the  leading  edge.  The  pressure  falls  from  this  value  at  the  Mach 
line  where  tan  0  =  1  to  0  at  the  Mach  line  where  tan  0  -  -1. 

For  a  semi-infinite  symmetrical  wedge  with  a  subsonic  leading 
edge  (fig.  2(b))  the  conical  pressure  field  is  given  by  the  following 

equation: 


P 


E.P. 


2  t/< 


n  ./tan2  A  — 1 


cosh 


— l 


tan  A—  tan  0 
|  1  -  tan  A  tan  0  | 


(2) 


The  pressure  coefficient  rises  from  zero  at  the  near  Mach  line  where 
tan  S  -  1  to  infinity  at  the  leading  edge  and  falls  again  to  zero  at 

the  far  Mach  line  where  tan  0  =  -1. 
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It  is  convenient  to  consider  the  foregoing  solutions  for  wedges 
as  sources  or  sinks  depending  on  whether  they  produce  positive  or 
negative  pressures.  Thus  when  a  flow  is  symmetrically  diverged  as 
"by  a  wedge,  the  changes  in  pressure  are  positive  and  can  be  considered 
the  result  of  a  pressure  source  coincident  with  the  leading  edge  of 
the  wedge.  When  a  flow  is  symmetrically  converged,  as  when  flowing 
over  the  maximum  thickness  position  of  a  diamond  profile  wing,  the 
changes  in  pressure  are  negative  and  can  be  considered  the  result  of 
a  pressure  sink  lying  along  the  midchord  line. 

The  pressure  fields  of  the  sources  and  sinks  which  shape  the 
wing  are  added  to  give  the  wing  pressure  field.  From  the  leading 
apex  two  sources  are  swept  back  along  the  leading  edges  to  infinity 
producing  wedge  boundaries  along  their  entire  extent  (fig.  3(a)).  To 
confine  these  wedge  boundaries  to  the  finite  length  of  the  leading 
edges  of  the  wing,  two  image  sinks  of  equal  and  opposite  magnitude 
are  superimposed  onto  the  portions  of  the  sources  lying  off  the  wing 
(fig.  3(b)),  thereby  producing  in  part  the  effect  of  a  tip.  Two 
semi-infinite  sinks  are  placed  along  the  midchord  line  at  right 
angles  to  the  flow,  and  two  image  sources  are  superimposed  on  those 
portions  of  the  sinks  lying  beyond  the  wing  tips  (fig.  3(c)).  The 
two  image  sources  combined  with  the  two  aforementioned  image  sinks 
fully  represent  the  effect  of  the  tips.  To  return  the  flow  at  the 
trailing  edge  to  its  original  direction,  two  sources  (fig.  3(d)) 
are  started  at  the  extremities  of  the  trailing  edges  and  are  extended 
back  along  the  trailing  edges  to  infinity  intersecting  at  the 
trailing  apex.  The  images  of  the  trailing-edge  sources  cannot  have 
any  contribution  to  the  wing  pressures  and  need  not  be  considered. 


Wing  Drag  Coefficient 

The  contribution  to  the  wing  pressure  drag  of  the  pressure 
field  of  each  source  and  sink  can  be  individually  determined,  and 
the  sum  of  the  individual  pressure  drags  will  then  be  the  total 
drag.  This  may  be  expressed  in  coefficient  form  as 

CD  -  Cpa  +  +  Cpc  +  cDd  +  cDe  (3) 

where  the  symbols  have  the  following  meanings: 

drag  coefficient  associated  with  the  leading-edge  sources 

CDb  drag  coefficient  associated  with  the  midchord— line  sinks 

Cp  drag  coefficient  associated  with  the  trailing-edge  sources 

c 
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CDd  drag  coefficient  associated  with  the  leading-edge  image  sinks 

CD  drag  coefficient  associated  with  the  midchord-line  image 

®  sources 

The  derivation  of  the  drag  coefficients  is  performed  for  a  wing 
of  type  8,  because  the  results  for  this  wing  type  can  he  reduced  to 
include  every  other  wing  type  as  a  special  case.  No  loss  of 
generality  is  incurred  by  restricting  the  analysis  to  M2  =  2  as 
will  be  subsequently  explained. 

The  actual  determinations  of  Cj^,  CpQ,  Cpd,  and  CpQ 

are  carried  out  in  Appendix  A.  The  expressions  for  these  drag  coeffi¬ 
cients  are  given  in  terms  of  t,  the  tangent  of  the  leading-edge 

sweep  angle,  and  §,  the  ratio  of  the  semi  span  to  the  root  chord.  For 
a  given  value  of  aspect  ratio  and  taper,  the  values  of  t  and  § 
are  determined  from  the  following  expressions: 


T  s  2(1-21 

A  (1+7) 


(4) 


RESULTS  AHD  DISCUSSION 


Design  Charts 

Because  the  equations  are  too  long  to  be  very  useful  for  direct 
computation  of  the  drag  coefficients,  detailed  calculations  have  been 
undertaken  to  determine  drag  charts  based  on  these  equations.  Although 
the  analysis  and  calculations  have  been  performed  for  M2  —  2,  the 
drag  charts  have  been  generalized,  as  subsequently  discussed,  to  be 
valid  for  an  supersonic  Mich  numbers. 

A  given  wing  operating  at  M2>2  or  M2<2  may  be  converted  to 
an  equivalent  wing  at  M£  *=  2  by  changing  its  lateral  dimensions  by 

the  factor*/  M 2  — 1  and  holding  the  thickness  ratio  constant.  The 

effective  aspect  ratio  of  the  equivalent  wing  is  A */M^— 1,  but  its 
taper  is  the  same  as  that  of  the  given  wing.  During  the  lateral 
expansion  or  contraction  the  relative  positions  of  the  Mach  lines 
and  the  wing  are  unchanged,  and  the  pressure  coefficients  computed 
for  the  equivalent  wing  at  if  =2  will  bear  a  point-to-point 
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correspondence  with  the  pressure  coefficients  for  the  given  wing  at 
M2>  2  or  M2<  2.  However,  the  pressure  coefficients  and  the  drag 
coefficients  for  the  equivalent  wing  must  he  divided  by  j/M2-!  to 
convert  them  hack  to  the  given  wing.  From  these  considerations  it 

follows  that  plotting  Cp  «/m2-1  against  A«/m2-1  for  any  taper 
will  give  a  chart  which  is  applicable  to  all  supersonic  Mach  numbers. 

The  drag  coefficients,  Cpa,  Cr^,  Cpc,  and  CpQ  have 

all  been  calculated  for  M2=2  and  then  have  been  added  to  give  Cp, 
the  wing  drag  coefficient.  Drag  charts  have  then  been  constructed 

to  give  — ~ —  vta2— 1  as  a  function  of  k>J  M2-!,  the  effective 

(|)2 

aspect  ratio,  for  various  tapers.  The  drag  chart  for 

0<  A  </M2~i£  1.8  is  given  in  figure  4(a)  and  the  drag  chart  for 

1.8<A  Vm2-1<  5.4  is  given  in  figure  4(b). 

An  examination  of  figure  4  reveals  several  points  of  interest. 

For  any  Mach  number  greater  than  unity,  the  drag  coefficient  decreases 
toward  zero  as  the  aspect  ratio  decreases  toward  zero.  Also,  As  would 
be  expected  the  drag  coefficients  approach  the  section  drag  coeffi— 

4(t/c)2 

cient  -  ■  w-  ■■  ...  as  the  aspect  ratio  Increases  toward  infinity.  The 

figure  also  shows  that  the  effect  of  taper  is  generally  greater  at 
small  effective  aspect  ratios  than  at  large  effective  aspect  ratios. 

Figure  4  also  reveals  the  following  information:  The  curves 
of  drag-coefficient  parameter  as  a  function  of  effective  aspect 
ratio  are  continuous  for  all  tapers.  However,  all  the  curvee  have 
discontinuities  in  slope  except  the  curve  for  X  =  1.  The 
discontinuities  in  slope  occur  when  the  leading  and  trailing  edges 
are  coincident  with  their  Mach  lines,  and  correspond  to  conditions 
given  by  line  III  of  figure  1.  Slight  decreases  in  effective  aspect 
ratios  below  those  for  this  condition  correspond  to  subsonic  leading 
and  trailing  edges  and  are  accompanied  by  a  rapid  decrease  in  the 
drag  coefficient,  particularly  for  small  tapers.  The  wing  having  the 
greatest  drag  coefficient  is  the  wing  of  zero  taper  and  effective 
aspect  ratio  2;  that  is,  the  wing  with  the  leading  edge  parallel 
to  the  Mach  lines.  This  wing  has  a  drag  coefficient  only  about 
15  percent  greater  than  the  section  drag  coefficient.  Below  an 
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effective  aspect  ratio  of  1.75,  wings  of  diamond  plan  form  have  the 
least  drag;  above  this  effective  aspect  ratio  rectangular  wings  have 
the  least  drag. 

As  an  example  in  the  use  of  these  charts  to  determine  the  wing 
drag  coefficient,  consider  a  wing  of  aspect  ratio  2,  taper  0.5,  and 
of  thickness  0.1  (same  as  semiwedge  angle  in  radians  for  thin  wings) 

flying  at  M  *»  1.3*  Since  A  1.662,  figure  4(a)  gives 

_ SB —  vfa2_l  =  4.23.  Thus  Op  =  (0.l)e  (1.204) (4.23)  =  0.051. 

(t)2 


It  can  he  noted  that  the  drag  curve  for  the  diamond  plan -form 
wing  may  he  used  for  wings  incorporating  a  flat  lower  surface  and 
an  isosceles  triangular  section  if  the  effective  aspect  ratio  is 
greater  than  two.  For  such  a  wing  with  no  tip  and  supersonic  leading 
edges,  the  upper-  and  lower-surface  pressures  are  independent  so  that 
the  drag  for  the  wing  of  isosceles  triangular  section  is  just  half 
the  drag  for  the  wing  of  the  same  plan  form,  the  section  of  which  is 
a  diamond  formed  hy  two  of  the  isosceles  triangles. 


Jfegnitude  of  Individual  Drag  Coefficients  for  Wings 

of  Rectangular  and  Diamond  Plan  Form 

Thft  general  formulas  for  the  individual  drag  coefficients  have 
been  reduced  to  their  particular  forms  for  rectangular  wings,  and 
the  results  are  given  in  Appendix  B.  The  ratios  of  the  individual 
drag  coefficients  to  thickness  ratio  squared  as  functions  of  aspect 
ratio  are  given  in  figure  5(a)  f  °r  M?=2.  The  leading— edge  sources 
produce  uniform  pressure  over  both  front  and  rear  surfaces  and  thus 
have  no  contribution  to  the  drag.  The  ridge  sinks  account  for  a 
drag  coefficient  equal  to  the  section  drag  coefficient  at 

all  Mach  numbers,  which  accounts  for  most  of  the  wing  drag. 

From  figure  5  (a)  it  may  also  be  seen  that,  above  an  aspect 
ratio  of  unity,  Cp^  the  drag  associated  with  the  leading-edge 

ineise  sinks  is  just  offset  by  Cpe  the  thrust  associated  with  the 

midchord-line  image  sources.  Thus  the  tip  does  not  affect  the  net 
wing  drag  although  it  does  affect  the  spanwise  distribution  of  the 
drag  near  the  tips.  Since  the  tip  contributes  no  drag  to  the  wing 
at  aspect  ratios  greater  than  1.0,  the  wing  drag  coefficient  is 
independent  of  aspect  ratios  in  this  range  and  must  be  equal  to  the 
section  drag  coefficient. 


MCA  TN  Wo.  1487 


9 


For  aspect  ratios  less  than  1.0  the  tip  affects  the  wing  drag 
in  a  manner  such  that  the  drag  is  reduced  because  the  negative  drag 
associated  vith  the  midchord— line  image  sources  more  them  offsets 
the  positive  drag  associated  with  the  leading-edge  image  sinks. 

When  the  aspect  ratio  approaches  zero,  the  thrust  of  the  images  is 
sufficient  to  overcome  the  drag  associated  with  the  midchord  line, 
so  that  the  drag  coefficient  also  approaches  zero. 

Also  in  Appendix  B  the  general  formulas  for  the  individual  drag 
coefficients  have  been  reduced  to  their  particular  form  for  wings  of 
diamond  plan  form.  These  results  have  been  used  as  the  basis  for 
figure  5(h)  wherein  the  values  of  the  ratio  of  the  individual  drag 
coefficients  to  the  thickness  ratio  squared  have  been  given  as  a 
function  of  aspect  ratio  for  V^-2. 

With  reference  to  figure  5(h),  it  can  he  seen  that  for  aspect 
ratios  greater  than  2.0  the  drag  associated  with  the  leading-edge 
sources  is  additive  to  that  associated  with  the  midchord  line  to 
make  the  wing  drag  coefficient  greater  than  the  section  drag  coeffi¬ 
cient.  For  aspect  ratios  less  than  2.0  the  wing  drag  coefficient 
falls  with  decreasing  aspect  ratio  as  the  result  of  several  effects. 
First,  the  trailing  edges  of  the  wing  are  now  subsonic  and  the 
positive  pressures  due  to  the  trailing-edge  sources  extend  onto  the 
wing  rear  surfaces  and  act  to  decrease  the  drag.  Second,  the  effect 
of  the  tip  is  to  decrease  the  drag  since  the  thrust  associated  with 
midchord  image  sources  is  greater  than  the  drag  associated  with  the 
leading-edge  image  sinks.  Finally,  the  drag  associated  with  the 
leading-edge  sources  decreases  because  the  sweep  of  the  sources 
increases  as  the  aspect  ratio  becomes  smaller. 

The  tip  drag  is  represented  by  the  difference  between  the 
drag  associated  with  the  leading-edge  image  sinks  and  the  midchord— 

line  image  sources.  For  a  rectangular  wing  with  A  1  >1 

the  drags  identically  cancel  each  other  so  that  the  tip  drag  is 
zero.  This  result  is  also  true  for  other  taper  ratios  as  long  as 
the  Mach  lines  of  the  images  intersect  the  trailing  edge  of  their 
own  half  wing.  It  has  also  been  found  in  other  calculations  not 
presented  in  this  report  that  for  tapered  and  un tapered  swept  wings 
with  a  diamond  profile  and  tips  cut  off  in  the  stream  direction  the 
drag  of  the  tip  is  zero,  provided  that  the  tip  of  one  half  wing 
does  not  affect  the  opposite  half  wing. 


Variation  of  Drag  Coefficient  with  Mach  Number 


Using  the  drag  charts  of  figure  4,  the  variation  with  Mach 


number  of 


has  been  determined  for  rectangular  wings 
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of  several  aspect  ratios  and  the  results  are  presented  in  figure  6(a). 
The  results  show  that  for  a  given  Mach  number  there  exists  a  certain 
aspect  ratio  above  which  all  rectangular  wings  have  drag  coeffi- 
dents  equal  to  the  section  drag  coefficient#  Thus  all  rectangular 
wings  of  aspect  ratio  greater  than  unity  have  the  section  drag 
coefficient  at  a  Mach  number  greater  than  1.41.  At  Mach  numbers 
near  unity  very  low-aspect-ratio  wings  show  some  reduction  in  drag 
compared  with  wings  of  higher  aspect  ratio. 

The  variation  with  Mach  number  Cj)/^  —  ^  has  also  been 

determined  for  wings  of  diamond  plan  form  of  several  aspect  ratios, 
and  the  results  are  presented  in  figure  6(b).  The  results  shew 
that  considerable  drag  reduction  is  to  be  realized  at  Mach  numbers 
near  unity  by  reducing  the  aspect  ratio. 


CONCLUSIONS 

The  following  conclusions  may  be  drawn  from  the  analysis  within 
the  limits  of  first-order  theory;  that  is,  insofar  as  the  disturbance 
velocities  are  small  and  the  effect  of  viscosity  may  be  neglected. 

1.  Below  an  effective  aspect  ratio  of  1.75,  wings  of  diamond  plan 
form  have  the  least  drag;  above  this  effective  aspect  ratio  rectangular 
wings  have  the  least  drag. 

2.  At  very  low  effective  aspect  ratios,  increasing  the  ratio 
of  tip  chord  to  root  chord  increases  the  drag  coefficient  markedly, 
but  at  large  effective  aspect  ratios  increasing  the  ratio  of  tip 
chord  to  root  chord  slightly  decreases  the  drag  coefficient. 

3.  Rectangular  wings  of  effective  aspect  ratio  greater  than 
unity  have  a  drag  coefficient  e^ual  to  the  section  drag  coeffi¬ 
cient. 

4.  For  wings  the  tip  Mach  lines  of  which  intersect  the 
trailing  edges  on  their  own  half  wing,  the  tip  effect  does  not 
change  the  net  drag  but  changes  only  the  distribution  of  the  drag. 

5.  With  decreasing  effective  aspect  ratio  for  a  given  taper 
there  is  a  sudden  decrease  in  drag  coefficient  as  the  leading  and 
trailing  edges  pass  through  the  Mach  cones. 

6.  At  low  effective  aspect  ratios  the  drag  coefficient 
decreases  with  effective  aspect  ratio. 

Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  For  Aeronautics, 

Moffett  Field,  Calif. 
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APPENDIX  A 

DETERMINATION  OF  INDIVIDUAL  DRAG  COEFFICIENTS 
Leading-edge  Sources 

Considering  first  the  leading-edge  sources,  figure  3(a)  shews 
the  five  areas  of  the  half  wing  over  which  their  pressure  fields 
must  he  Integrated  to  determine  the  drag  coefficient  Cp^.  Since 
there  are  two  half  wings  and  two  sides  to  each  half  wing,  and  since 
the  front  surface  is  inclined  at  +t/c  radians  to  the  flow  and  the 
rear  surface  Is  Inclined  at  — t/c  radians  to  the  flow 


Considering  both  leading-edge  sources,  the  pressure  coefficient  is 
given  by  equation  (2)  as  follcvs: 

+  cosh"1  I— |)  (A2) 

1+T0  [  I  1-T0  IJ 


R.P. 

* 


_2t/o 

/ 


-fooeh" 

T2— 1  ^ 


The  differential  areas  taken  as  small  triangles  from  the  apex  or  as 
differences  of  the  triangles  are  as  follows: 


dSi  «*  — -  d a 

8 


ds2  =  —  da 
8 


dS 


3"  25? 


(A3) 


sf. 
2 a1 


c* 


dSA  =  ^ —  da  —  —  da 

+  -  .2  3 


dS 


da 


da 


2  (l+o  t)2  8 
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Substituting  from  equations  (A2)  and  (^3)  into  aquation  (Al), 
putting  in  the  limits  shown  in  figure  3(a) j  and  simplifying 
gires  the  following  result: 


28 


c-e  t 


Jjjosh  1 


T+q 

1+Tq 


+  cosh  1 


r-q 

1-to 


Carrying  out  th»  Integration,  collating  einilar  terte,  mking 
the  following  eutatltutlon  for  the  wing  area  S  -  S  so  (l-e  /o) 


gives  the  drag  coefficient. 
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Cd&  *=  E.P. 


(2b/c) 


+  ft+jW«>g  ^ / °  1  _  fr-e°V°r.  c08h- 

T  L  1+28  t/c  J  T 


- 

L  l-2s  t/c  j 


+  j.l-7Tg)+4(T2-l)sT/c  C08h-1[T_B(T2-1)/C] 

T(l— T2) 

+  (l-ggT/c)2  cosh-1  r-^T2+l)/C  -  — 2t2~1-^  cosh-1  T 

T  1— 2bt/c  •*  t(  T2— l) 


2 


«A  2_1 


2  </  (1-bt/c)2-(b/c)2 


4s2t 

c2 


(A5) 


Midchord— line  Sinks 


A  simple  expression  can  be  derived  for  the  drag  coefficient 
associated  with  the  pair  of  sinks  along  the  mldchord  line.  From 
equation  (1)  the  pressure  drop  for  the  flow  deflected  through  an 


angle  -2 


(0 


by  the  sinks  is 

4 


*b 


(I) 


4  if ;  r  1 

- - \./.  COS  x(-o)  +  cos  1  (o) 

ir 


Since  the  pressure  coefficient  acts  uniformly  over  the  upper  and 
lower  rear  surfaces  to  cause  drag. 


<*>  -  *  (IT 


(A7) 


Ik 
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Trailing— edge  Sources 

The  sources  along  the  trailing  edges  (fig.  3(d))  act  always 
to  decrease  the  drag  as  their  positive  pressure  field  may  act  only 
on  the  rear  surface  of  the  wing.  The  drag  coefficient  is  given  by 
the  following  equation: 


pc(|)asi 


The  pressure  coefficient  for  one  source  from  equation  (2)  is 


Pc  =  R.P. 


2t/c  . 

— - —  cosh 


nJK 1 


and  the  differential  area  is 


dS  =  2s*  T' 


(1+T0) 


(A10) 


Substituting  from  equations  (A9)  and  (A10)  into  equation  (a8), 
carrying  out  the  integration,  and  simplifying  gives  the  drag  coeffi¬ 
cient 


.  /t\  /s\ 

+  l TS  )  V5  )  r 

Cn  =  -E.P.  — VZ—  v  z - 

°  It  V*-2— 1  (1-flT/c) 


.  _x  |T +1 

cosh  1  j-gr- 


Leading-edge  Image  Sinks 

The  image  sinks  associated  with  the  leading— edge  sources  produce 
pressures  on  the  wing  which  tend  to  increase  the  wing  drag.  Referring 
to  figure  3(b),  the  drag  coefficient  can  be  written  as  follows: 

s  “  X Pd  (O43*  -  *  i>  (O'*  -  *.4. Pt  (0433 


-kf  pd(|)^  +  “£pd(|)^ 


(All) 
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The  differential  areas  are  given  by  the  following  equations: 


*2c2 

dSi  =  da 


A2c2  do  X2c2 

dS2  = - da 

2  (1-to)2  8 

*2c2 


dS  =  — 


do 


3  2  (1+To)2  8 


da 


(A12) 


dS4  =  2s' 

a * 


da  ■\scs 
_ 


do 


dS5  .  2s2  & 
5  a2 


The  pressure  coefficient  can  be  obtained  by  an  application  of 
equation  (2). 

Pd  =  -  B.P.  cosh"1 

K  V  T2— 1  I  1+TO 


(A13) 


Substituting  from  equations  (A12)  and  (A13)  hack  into  equation  (All), 
putting  in  the  limits  shown  in  figure  3(b),  and  simplifying  gives 
the  drag  coefficient. 


S  CD,  =  R.P, 


(t) 


2 

C2 


¥ 


As 


^  /c— esT  cosh 


— l 


T+O 

1+TO 


da 


2  P  ^ 

+  2  A  /  c_eT  cosh-1 


T+a 

1+TO 


da 


(1— to) 


2 


$ 


28 


+  2  /C—2  8T  cosh~i 


3 

C— -ST 


T+a 


1+Ta 


do 


(L  +tctY 


8s2/"  48  >.-1 

— —  /  C-2ST  cosh  1 

r+a  1 

da 

__  8sfr 1 

co=h"i  f  T+<?  N 

\  da  ~ 

J  2S 

C“68T 

1+Ta| 

a2 

~  C 2  7 

C—28T 

\1+Tav 

'a2  . 

(A14) 
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Carrying  out  the  integration,  collecting  similar  terms,  and 
substituting  S  =  2sc(1-st/c)  gives  the  drag  coefficient. 


rtN^-l(l-ST/c) 


(1+2st/c_)2  ooeh-1 


t— 2b(t2— 2)/c 
-  l+2sr/c 


(1-2st/c)2 


4s/c  (l— 2a  r/c)g 

1-2st/c  t 


cosh-1 


T 


+  C (31^—1)  +4st(1— T6)  /c j  coah_1  [t_3  (t2_d  /c  ] 

T  (re— 1) 


+  (l-28  T/°j_  Cosh  1 

T 


t-b  (t2+1)  /c 
1-28t/c 


_  2  [t-4s  (ts-1)  /cj  cQBh_1  jT  _2g  (t2_i}  /c  ]  +  2 
(t2— 1) 


4/  (1— 2st/c)2— (2s/c)g 


Z«/  (1— st/c)2— (a/c)g  +  16s2 

c2 


8s2 

c2 


cosh 


— l 


(A15) 
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Midchord-line  Image  Sources 

The  image  souroes  associated  with  the  midchord— line  produce  positive 
pressures  on  the  rear  surfaces  of  the  wing  and  tend  to  reduce  the  drag. 
Referring  to  figure  3(c),  the  drag  coefficient  can  be  written  as  follows: 


s  £*,  =  -  k  p.  (l)  dSx  -  k  p.  (l)  -  if  P.  (l)  «Ss  (*16) 


The  differential  areas  can  be  written  as  follows: 


dSi 


*2c2  do 
~8  (l»ro)  2 


dS2 


c£ 

8 


(1+2st/c)2 


do 

(l+TO)2 


(A17) 


dS3  =  2s2  $2- 
o2 

The  pressure  coefficient  may  be  obtained  by  an  application  of  equation  (1) . 

Pe  =  fo/g.  cos-i  (a)  (Al8) 

It 

Substituting  from  equations  (A17)  and  (Al8)  into  equation  (Al6),  putting  in 
the  limits  shown  in  figure  3(c),  and  simplifying  gives  the  drag  coeffi¬ 
cient. 


SCd6 


(a) 


do 

(I-to)2 


(a) 


do 


(1+to) 


i6sg  r1 
,2jL± i. 


c-asr 


cos 


(o) 


do 

o2 


(A19) 
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Carrying  out  the  integration,  collecting  similar  terms,  and  making  the 
substitution,  S  =  28c(1-st/c),  gives  the  drag  coefficient. 
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APPENDIX  B 

REDUCTION  OF  GENERAL  EQUATIONS  FOR  WINGS  OF 
RECTANGULAR  AND  DIAMOND  PLAN  FORM 

The  general  formulas  for  the  individual  drag  coefficient  can  be 
considerably  simplified  for  vings  of  rectangular  or  diamond  plan  foim. 
For  the  rectangular  plan  form  the  formulas  are  simplified  by  determin¬ 
ing  the  limiting  values  of  the  drag  coefficients  as  t,  the  tangent 
of  the  leading— edge  sweepbeck  angle,  approaches  zero  since  the  formulas 
are  indeterminate  for  t=  0.  The  actual  limiting  process  is  lengthy 
and  only  the  results  are  given  here. 


Summing  the  component  drag  coefficients  and  making  the  substitution 
A  =  2s/c  gives  the  wing  drag  coefficient  in  terms  of  the  aspect  ratio. 
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CD 


+  A  cosh 


1 

A 


—  2  cos 


-X 


V 


_  i  (—  tL - 4a  cosh  1  —  +  4  cos  1  2 a'') 

n  V  A  2A  / 


(B2) 


Inspection  of  equation  (B2)  shows  that  the  first  term  is  real  for  all 
aspect  ratios,  that  the  next  three  terms  are  real  only  for  aspect  ratios 
less  than  unity,  and  that  the  last  three  terms  are  real  only  for  aspect 
ratios  less  than  one-half . 

In  simplifying  the  general  formulas  for  wings  of  diamond  plan  form 
it  is  sufficient  to  substitute  the  condition  for  zero  taper  r  =  c/2s, 
since  the  formulas  are  determinate  for  this  case.  The  following  results 
have  been  obtained. 


CDa  =  E.P. 


(sX'i 

r2(2rs-l) 

•J^-1  ' 

It(1-t2)  _ 

cosh  1  t  —  cosh" 


NACA  TN  No.  1487 


21 


Summing  the  component  drag  coefficients  and  making  the  substitution 
A  =  2/t  gives  the  following  formula  for  the  drag  coefficient  of  a 
diamond  plan— form  wing: 


CD  =  E.P. 


f  a 

It  '  */4-A2 


1  +  -  *J  h-Az  sin”1 
2 


,  8-Aa^ 
\k=^J 


2  cosh- 


4+A£ 
V  4A 


) 


-  cosh 


(B4) 


When  the  leading  edges  of  the  wing  are  swept  behind  the  Mach  cone, 
r  >1  and  A<  2  so  that  all  terms  of  equations  (B3)  and  (B4)  are  real 
However,  for  leading  edges  in  front  of  the  Mach  cone  where  t  <  1  and 
A>  2,  use  must  be  made  of  the  following  relationship  to  1 ind  the  real 

part  of  equations  (B3)  and  (B4) : 


cos  1  x  =  cosh  1  x  (bj) 

•J  1_-T2  */t2— 1 

The  resulting  equations  for  the  component  drag  coefficients  and  wing 
drag  coefficient  are  as  follows: 


2  .  _  £ c08— l  T 

yil^2  Td-^) 


CD 


4 

16  %  , 

m  .  . . 

it  v  A2— 4 


-  a2— 4  — 

.  2 


A 

Va2-i> 


('S')  1 

2  (A2-4)  V'A/  j 


) 

i 

.i 

I 

» 

>  (B6) 


i 

! 

j 

! 

i 

\ 

\ 

J 


(BT) 
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Fig.  1 


Taper,  A 

Figure  /.-  Subdivision  of  wings  into 
eight  types  for  M*  -  2 


Fig.  2  a,  b 
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Leading  edge  (9  + A  s  -Jr  ) 
/  Mach  line  (  tan  9  -  /  ) 

A 


M*=  2 


-A.  (  45  deg. 


\  Mach  line  ( ton  9  =  -  / ) 

~2ib 
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Section  A  A 


(a)  Supersonic  leading  edge 


,  Mach  line  ( tan  9  -  I  ) 

/ 


(b)  Subsonic  leading  edge 
Figure  2.-  Symmetrical  semi~ infinite  wedges  of 
Jones  basic  solution 
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Fig.  3  a 


Figure  3-  Sources,  sinks ,  integration  areas,  and  limits 
for  determining  wing  drag  coefficient 
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(b)  High  aspect  ratio  range 
Figure  4-  Concluded 


Ratio  of  drag  coefficient  to  thickness  ratio  squared 
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Fig.  5  a 


Figure  5-  Ratio  of  Individual  drag  coefficients  to  thickness 
ratio  squared  as  a  function  of  aspect  ratio  at  M2  -  2 


Ratio  of  drag  coefficient  to  thickness  ratio  squared 


Fig.  5  b 
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Aspect  ratio ,  A 
(b)  Diamond  ptanform 
Figure  5r  Concluded 


Ratio  of  wing  drog  coefficient  to  thickness  ratio  squared, 
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Two-dimensional  (Ackeret) 
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Mach  number ,  M 
(a)  Rectangular  pianform 
Figure  6—  Ratio  of  wing  drag  coefficients  to  thickness 
ratio  squared  for  several  aspect  ratios  as  a  function 
of  Mach  number 


Ratio  of  wing  drag  coefficient  to  thickness  ratio  squared , 


Fig.  6  b 
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Mach  number,  M 


(b)  Diamond  pianform 
Figure  6-  Concluded 


